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A simple set of equations is derived for predicting the temperature distribution and radiative energy flow in 
a semitransparent layer consisting of an arbitrary number of laminated sublayers that absorb, emit, and scatter 
radiation. Each sublayer can have a different refractive index and optical thickness. The plane composite region 
is heated on each exterior side by a different amount of incident radiation. The results are for the limiting case 
where heat conduction within the layers is very small relative to radiative transfer, and is neglected. The 
interfaces are assumed diffuse, and all interface reflections are included in the analysis. The thermal behavior 
is readily calculated from the analytical expressions that are obtained. By using many sublayers, the analytical 
expressions provide the temperature distribution and heat flow for a diffusing medium with a continuously 
varying refractive index, including internal reflection effects caused by refractive index gradients. Temperature 
and heat flux results are given to show the effect of variations in refractive index and optical thickness through 
the multilayer laminate. 
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Nomenclature 

absorption coefficient of yth sublayer, m’ 1 
total thickness of the composite plane 
layer, m 

thickness of yth sublayer in composite, m 
exponential integral functions, 

= LV m - 2 exp (— j c//a) d n 
function of refractive index defined in 
Eq. (19a) 

index of refraction of yth sublayer 
heat flux, W/m 2 

externally incident radiative flux, W/m 2 
dimensionless radiative heat flux, 

qAq°n ■ q°2) 

quantity pj/(l - pj) 
absolute temperature, K 
dimensionless temperature function, 

(ctT 4 - q 0 r2 )l(q°, ~ q° 2 ) 
dimensionless local coordinate in yth 
sublayer, XjIDj 

coordinate normal to boundaries of 
composite plane region, m; x Jy local 
coordinate in yth sublayer, m 
optical thickness of entire layer 

+ k 2 + *•* + Kj\ K Dj optical thickness of 
the yth sublayer (a } + (t S] )Dj 
optical depth in a sublayer ( a } + v SJ )Xj ; 
extinction coefficient in the complex index of 
refraction 

radiation wavelength in vacuum 
reflectivity of interface for internally and 
externally incident radiation 
Stefan-Boltzmann constant, W/(m 2 -K 4 ) 
scattering coefficient in yth sublayer, m" 1 


Received Sept. 4, 1992; revision received Nov. 23, 1992; accepted 
for publication Nov. 24, 1992. Copyright © 1992 by the American 
Institute of Aeronautics and Astronautics, Inc. No copyright is as- 
serted in the United States under Title 17, U.S. Code. The U.S. 
Government has a royalty-free license to exercise all rights under the 
copyright claimed herein for Governmental purposes. All other rights 
are reserved by the copyright owner. 

^Research Scientist, Lewis Research Academy. Fellow AIAA. 
tResearch Scientist, Heat Transfer Branch. 


r', t° = transmissivity of interface for internally and 
externally incident radiation 
= dimensionless temperature function, Eq. (3) 
W = dimensionless radiative heat flux function, 

Eq. (4) 


Subscripts 

H , L 
/, o 
J 

i 


= first (Xj = 0) and second (x, = D } ) internal 
interfaces of a sublayer within the composite 
= higher and lower values 
= incoming and outgoing radiation 
= total number of sublayers in the composite 
= yth sublayer in the multilayer composite, 

1 < j < J 

= radiative quantity 


Superscripts 

i = inside surface of an interface 

o - outside surface of an interface 

Introduction 

T HE use of ceramic parts or coatings is of interest for 
high- temperature applications. Some parts may have 
reinforcing layers or may be laminated, so it is necessary to 
consider heat transfer in composite materials. Surrounding 
temperatures are usually high, so there can be appreciable 
heating by radiation. Since some ceramics are semitranspar- 
ent, their temperature distributions are influenced by internal 
radiative heat flow. The refractive indices of the sublayers in 
a composite can have a considerable effect on the temperature 
distribution and radiative heat flow. Surface reflections de- 
pend on the refractive indices on both sides of the interface, 
so the refractive indices influence both the external energy 
transmitted into the interior of a material and the amount of 
internal energy reflected back into the interior. Within a ma- 
terial, emission depends on its refractive index squared; hence, 
internal emission can be many times that for a blackbody 
radiating into a vacuum. Since internal radiation exiting through 
an interface into a vacuum cannot exceed that of a blackbody, 
there can be considerable reflection at the internal surface of 
an interface, partially by total internal reflection. If there is 
a gradient of refractive index in the material, it will cause 
some of the radiation to be refracted backward when initially 
directed into a region with a smaller refractive index. These 
effects distribute internal energy within the material and make 
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the temperature distribution more uniform than in a region 
with a refractive index close to unity. This article will provide 
some insight into these effects by examining the limiting case 
where heat transfer is only by radiation. 

In Siegel and Spuckler 1 * 2 it was shown that for radiative 
equilibrium (no conduction or convection) in a single- or two- 
layer gray medium with diffuse interfaces, the temperature 
distribution and radiative heat flux for any refractive index 
can be readily obtained from available results for a single 
layer with a refractive index of one. This article uses these 
ideas extended to a multilayer laminate to determine radiative 
transfer characteristics in a layer with a continuous variation 
of refractive index. Analytical expressions that are readily 
evaluated are obtained for the temperature distribution and 
radiative heat flow. Each sublayer emits, absorbs, and iso- 
tropically scatters radiation. For simplicity, the external me- 
dium on each side of the entire laminate has a refractive index 
of one. 

The outer surfaces of the laminate, and all interfaces be- 
tween adjacent layers, are assumed diffuse. This is probably 
a reasonable approximation for ceramics that have not been 
polished and are bonded together. When transmitted radia- 
tion, or radiation emitted from the interior, reaches the inner 
surface of an interface, it is assumed diffuse as a result of 
scattering within the medium. If the index of refraction of the 
material is greater than that of the surrounding medium, some 
of this diffuse energy is in angular directions for which there 
is total internal reflection. These reflections provide a trap- 
ping effect that retains energy within the layer with the larger 
refractive index, and tends to equalize the local temperature 
distribution in this layer. 

There have been many studies of radiative transfer in semi- 
transparent layers. However, very little has been done to 
examine the effects of a variable refractive index, especially 
with regard to computing internal temperature distributions 
when there is internal emission and scattering. To predict heat 
treating and cooling of glass plates, Gardon 3 developed an 
analysis for temperature distributions in absorbing-emitting 
layers, including uniform refractive index effects. The inter- 
faces were optically smooth, and specular reflections at the 
interfaces were computed from the Fresnel reflection rela- 
tions. A similar application, Fowle et al., 4 predicted heating 
in a window of a re-entry vehicle. Some recent papers 5-8 
further examined the effects of Fresnel boundary reflections 
and a nonunity uniform refractive index. Many analyses of 
both steady and transient heat transfer to single or multiple 
plane layers, such as Amlin and Korpela 9 and Tarshis et al., 10 
have used diffuse conditions at the interfaces as in the present 
study. Thomas 11 set up a solution procedure to include a 
ceramic interface that is partially specular and diffuse. Ritchie 
and Window 12 and Snail 13 analyzed the reflection behavior of 
thin films with a graded refractive index. The results are con- 
cerned with the effect of thin films on the absorption of solar 
radiation and do not consider their internal temperature dis- 
tribution. 

Analysis 

A laminated plane layer consists of an arbitrary number of 
] sublayers of different materials with thicknesses D,, D 2 , 

. . . , Djy . . . , Dj, as in Fig. la. Each sublayer absorbs, 
emits, and isotropically scatters radiation. The limiting case 
is examined here where energy transfer within the laminated 
region is dominated by radiation, so heat conduction is ne- 
glected (in the literature this limit is referred to as “radiative 
equilibrium”). An analytical solution is obtained which can 
provide a helpful basic limiting case for numerical studies in 
which conduction is included. Each sublayer has a constant 
rij\ it is the effect of the different refractive indices and the 
different optical thicknesses of the sublayers that is investi- 
gated here. The materials are assumed to provide significant 
scattering, so the energy fluxes are assumed diffuse at the 



Fig. 1 Laminated multilayer geometry: a) coordinate systems in lay- 
ers, and nomenclature designating interfaces and their heat fluxes; b) 
detailed nomenclature for the jth layer. 

interfaces between sublayers and between the outer sublayers 
and the surrounding air or vacuum. For simplicity, the re- 
fractive index of the surroundings is one. 

As shown in Fig. la the laminated region is subjected, to 
diffuse q° rX and q° r2 from the external surroundings onto the 
two outer boundaries x x = 0 and x } — D } \ for convenience 
q°r\> q°i' Within each sublayer, the analysis will use outgoing 
and incoming fluxes, q o j and q i h that are moving away from 
or toward each interior interface of each sublayer. Since scat- 
tering is included, the local optical depth within each sublayer 
is related to its individual x f coordinate by k ; = {a } -I- cr Sf )x i 

= («; + aJ D i X j = «Di X i- 

The temperature distribution in each sublayer is governed 
by an integral equation (Siegel and Howell 14 ) with index of 
refraction factors included to give for the / th sublayer 

n)<rT](X r K Dj ) = | {q ojf E 2 (K Di Xj) + q OJ , s E 2 [K Dj ( 1 

- *,.)]} + | k Di l' a77(*;)£,(K 0/ |*/ - *;|) dx; 

( 1 ) 

where q o jJ and q OJrS are the outgoing internal radiative fluxes 
from the first and second interfaces of the j th layer as shown 
in Fig. 1. The energy equation d^ r /djc = 0 shows that q r being 
transferred through the entire laminated layer is constant for 
the present conditions of radiative equilibrium (negligible in- 
ternal heat conduction). The relation of q r to the temperature 
distribution in any of the sublayers is obtained by using the 
expression for q r within a plane layer as given in Ref. 14. This 
is evaluated at x = 0 to give for the j th sublayer 

Qr{ K Dj) = Qo,j,f ~~ ^io,J,s^3( K Dj) 

- 2 njK Di | 0 ‘ ctTHXjWkojXj) dX, (2) 
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The following dimensionless groups are now defined for each 
sublayer: 


(lsjs;) (3) 


Vj(l<Dj) = 


<?r(«D,) 


(1 s/sj) 


Qo,j,f QoJj 

In terms of these quantities Eqs. (1) and (2) become 


(4) 


Qr Qoj,f Qi,j,f 

Qo,j,f ~ Qij-l,s T j,f + 

The q i jtf is eliminated to yield 

Ro.i.f = [1/(1 “ P)Mq,.i-x. s -rlf - q r p), f ) (10a) 

At the second internal interface of the j - 1 layer, j - 1, s , 
the radiative flux is written in terms of incoming and outgoing 
fluxes as 


Xdj) = \ E 2 (K Di Xj) + | k di £ k d ,) 

x E,(k D) \Xj - X*\) dX* ( 1 S|< /) ( 5 ) 


Qr Qo,j-\,s 

This is used to eliminate q iyi -i, s fr° m Eq. (10a) which gives 
the following relation between the outgoing fluxes and the 
radiant flux q r being transferred: 


%(x Di ) = 1 - 2 k D/ 


, f *,«. 


k di )E 2 (k D jXj) dXj 


= [1/(1 - Pj,f)](<loJ-l.,Tj.f ~ <lrP‘j.f + qrTf.f) (10b) 


(i s/s;) (6) 

Equations (5) and (6) are valid for each of the J sublayers. 
It is important to observe that <!>, and 'Fy are not functions of 
rtj. Hence, to obtain and for all n i > 1, it is necessary 
to solve Eq. (5) numerically, only once, for each value of the 
parameter k D/ and use each result to evaluate Eq. (6). This 
has already been done, and tabulated results for <J>(AV k d ) 
and 'J'(ko) from numerical solutions are in references such 
as Heaslet and Warming 15 and Siegel and Spuckler. 1 ’ 2 

Although the numerical solutions for k d ) and 
are available, the solutions required here for q r through the 
composite layer, and for the Tj{X } ) in each of the sublayers, 
have not yet been obtained. This is because Eqs. (3) and (4) 
contain unknown outgoing fluxes q 0 jj > an d at the first 
and second interfaces of each sublayer. These boundary fluxes 
for each of the j sublayers must be expressed in terms of known 
quantities. This is accomplished by looking in detail at the 
outer boundary and interior interface conditions. 

At the two outer boundaries of the composite laminated 
layer, the internal fluxes leaving the internal interfaces of the 
first and / th layers are related to the transmission of external 
flux from the surroundings and the reflection of incoming 
internal flux by 

<7o. 1 ./ = q°r\T°, f + Qi.\.fP‘\.f (7a) 

q 0 .j. s = qfrls + <iu.spj.s (7b) 

The superscripts o and i designate properties for radiation 
incident on the outside or inside of an interface. Thus r° j f is 
on the outside of the first interface of the j th sublayer as 
shown in Fig. lb. At the inside surfaces of the two outer 
boundaries, the following relations exist between the radiative 
flux being transferred and the outgoing and incoming fluxes: 

<7r = q a ,u - qiA.f (8a) 

q r = -qo,j,s + qu.s (8b) 

Equations (7a) and (8a) are combined to eliminate q i l f and 
similarly Eqs. (7b) and (8b) are used to eliminate q i J s . This 
yields 

q 0 .i, f = [i/(i - pi, f)](q°n Ti, f - q r p\, f ) (9a) 

q 0 .i., = [1/(1 - PjJKq^Th + qrPj, s ) (9b) 

For the j th layer, at the first internal interface /, /, the two 
equations that involve reflected and transmitted fluxes are 


As a result of total reflections at interfaces when radiation 
passes into a medium with a lower refractive index, there are 
the following relations at the first and last internal interfaces 
of the entire layer and at ^he first internal surface of the j th 


sublayer 2 * 16 : 

T?,/ = (1 - Pi./)«1 

(Ha) 


t°j, = (1 - P‘j,s)nj 

(lib) 


T lf = (! - Pj.f)( n i ln i- 1) 2 

(11c) 

Equation (11) 
(10b) to yield 

is used to eliminate the r from Eqs. 

(9) and 


q 0 .\,f = 

(12a) 


q 0 ,j,s = n jq°2 Rj. s q r 

(12b) 

q 0 .i, f = ( n ji n j-i) 2 q 0 ,j-i,s - qARj.f - («>/«> _ ,) 2 ] 

(12c) 


where Rj = p* /( 1 — p' ). Equation (4) is now rearranged and 
used in Eqs. (12) to eliminate the q o j s fluxes at the second 
interface of each sublayer. This yields a set of simultaneous 
equations in terms of the q Q ,j,f fluxes at the first interfaces, 
the radiative flux q r and the function. These equations are 
arranged into the following array: 
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The equations in (13) are now added, which eliminates all of 
the q a terms; this yields the equation for the radiative flux q T , 
as 



reflected energy integrated over all incident directions for a 
diffuse incident distribution gives the relation for p'(/i) 14 


P'(n) s F(n) = | + 


(3/1 + 1)(/1 - 1) 
6 (n + l) 2 


n 2 (n 2 - l) 2 / /t — l \ 2 n 3 (n 2 + 2/i - 1) 

+ (n 2 + l) 3 + 1/ (« 2 + 1)(« 4 - 1) 


8n 4 (n 4 + 1) 

+ ( n 2 + 1)(« 4 - l ) 2 


A(«) 


(19a) 


The temperature distributions in the individual sublayers 
are now obtained. Using the set of Eqs. (13), the first j equa- 
tions are added to obtain 



(15) 


This is for diffuse radiation incident on a material with a higher 
refractive index, where n = n H ln L ( n H and n L are the “higher” 
and “lower” n values). When diffuse radiation is traveling in 
the reverse direction from a higher to a lower n value material, 
p' is given by 16 

P'(h) = 1 - (1/« 2 )[1 - f(«)] ( n = n H /n L ) (19b) 


The desired temperature distribution 7}(A)) is now written in 
terms of a dimensionless function t*(Xj) as 

7-4 _ Z 'll Z Zl 2 = zIzEl ~ zhh Sl 
' q°r i - q ° 2 q, «/(<?? i - qh) 

_ <?. / nj<rTf - + q - njq% \ 

q°n - q°r 2 \ njq r njq r ) 


The ratio of Eqs. (3) and (4) is substituted into the first term 
in parentheses of Eq. (16) to obtain 


TM _ 
1 J 


q°rl - q°r 2 


(*L 

W*< 


qpj.s 


n]q°r2\ 


(17) 


The negative of Eq. (15) is now inserted into the parentheses 
on the right side of Eq. (17). After rearrangement and using 
Eq. (4), the result for the temperature distribution (in terms 
of the dimensionless function ?}) of the ;th sublayer is 



<1</</) <18l) 


The first and last sublayers are special cases, and by similar 
manipulations their dimensionless temperature function dis- 
tributions are obtained from 


T 4 = 1 




l - <J>A 

i ) 


(18b) 


The reflectivity relations in Eq. (19) are derived for surfaces 
of dielectrics in the limit of infinite electrical resistivity, and 
for this condition they do not attenuate radiation internally. 
As discussed by Cox, 17 in the spectral regions where ceramic 
materials are reasonably transparent to radiation, so that in- 
ternal radiation will affect the temperature distribution, the 
extinction coefficient k (not to be confused with its optical 
thickness) in the complex index of refraction ( n - i#c) is 
usually not large enough to significantly affect the surface 
reflectivity. Hence, Eq. (19) for nonattenuating dielectrics 
often provides reasonable reflectivity results for attenuating 
dielectrics. The absorption coefficient a in a material, is re- 
lated to its extinction coefficient k by a = 47tk/A 0 . Since 
wavelengths for thermal radiation are in the micrometer range, 
only a small value of k is required to yield a large value for 
a. If k is large enough to influence the interface reflectivity 
relations, a = Attk/\ 0 will be so large that the radiating layer 
is essentially opaque, unless its thickness is much smaller than 
the ceramic layers considered here. 

The <!> and3^ are the building blocks to evaluate the present 
solution for a laminated composite layer that includes refrac- 
tive index effects and interface reflections. To evaluate results 
from Eqs. (14) and (18), the exact functions Opf, k d ) and 
¥(#c d ) can be used as obtained from numerical solutions of 
Eqs. (1) and (2); these functions are available in the literature 
as referenced earlier. Exact functions from the numerical in- 
tegration are tabulated in Siegel and Spuckler 1 for optical 
thicknesses, 0.1 ^ k d ^ 100. A convenient alternative to use 
for the functions ^(A^, k d ) and 'V(k d ) are the diffusion so- 
lution results which have the advantage of being simple an- 
alytical expressions. The use of these approximate functions 
to build up solutions was found to give accurate results for a 
two-layer composite. 2 The diffusion functions for a plane layer 
are given by 14 


Tj = 



(** + *L-) 

\nj nj%) 


(18c) 


The temperature distributions in all the layers can therefore 
be evaluated from Eqs. (18) after the dimensionless radiative 
flux q r /(q° r \ ~ q° r2 ) has been obtained from Eq. (14). 

To* use these relations, the internal interface reflectivity 
p *. values are needed to evaluate the Rj for various refractive 
indices of the sublayers. In the absence of better information, 
the following relations are used. The externally incident ra- 
diation is diffuse, and as a result of internal scattering and 
diffuse emission it is assumed that the internal radiation is 
also diffuse at the interfaces. Although the interfaces are not 
optically smooth, it is assumed that each bit of roughness acts 
as a smooth facet so that the reflectivity can be obtained from 
the Fresnel interface relations for a dielectric medium. The 


>P(k d ) = [1 !{\k d + 1)] (20a) 

0(*, k d ) = V(K D )[iK D ( 1 - X) + i] (20b) 

Very good comparisons were obtained for the present re- 
sults evaluated by using both the exact and the diffusion func- 
tions for d>(Af, k d ) and ^(#c D ). For this reason, and because 
tabular results for the exact solution are limited, the diffusion 
solution is recommended for simplicity and sufficient accuracy 
for heat flux and temperature distribution calculations using 
the present analytical equations. When used in Eqs. (14) and 
(18), Eqs. (20) provide rapid predictions of heat flux and 
temperature distributions for any refractive indices and op- 
tical thicknesses of the individual sublayers in the laminate. 
For the graphical results shown here, ^(A", k d ) and ^(k d ) 
were obtained from Eqs. (20), but the exact functions in Ref. 
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1 can also be used to yield similar results as shown on one of 
the figures. 

Results and Discussion 

Effect of Variation of Refractive Index 

Figure 2 shows the effect on the dimensionless temperature 
function of having the refractive index increase across the 
thickness of a plane layer. In Fig. 2a, the refractive index 
increases from one to three in equal increments within the 
five sublayers comprising the total thickness; in Fig. 2b the 
number of sublayers is increased to 21 to provide a more 
continuous variation of n over the same range of n. A large 
n increases the local internal radiation emission which is pro- 
portional to n 2 for a gray material. A larger refractive index 
in a material relative to its surroundings also increases the 
amount of internal reflection at its bounding interfaces. This 
is the result of total internal reflection when radiation is di- 
rected out of a material with a higher refractive index. These 
effects of refractive index act to distribute energy across a 
layer by virtue of emission and internal reflections. As a re- 
sult, in Fig. 2 the temperature profiles become more uniform 
as the n value increases across the layer thickness. 


so that the temperature in each sublayer is almost uniform, 
as is characteristic of an optically thin region. There are tem- 
perature jumps between layers. These temperature discon- 
tinuities occur at interfaces that separate regions of different 
refractive indices when heat transfer is only by radiation (the 
limiting condition analyzed here). The presence of heat con- 
duction would eliminate the temperature discontinuities at 
the interfaces, but there can be large local temperature gra- 
dients near the interfaces if conduction is small, relative to 
radiation. As shown in Fig. 2, the discontinuities decrease as 
the optical thickness increases. For k d = 100, the temperature 
profile in Fig. 2b is essentially continuous. 

As optical thickness increases, the resistance to radiative 
heat flow increases, so that the temperature gradients increase 
throughout the entire layer and within each sublayer. For k d 
= 0.1 the layer is optically thin, and the change in the di- 
mensionless temperature function across the entire layer is 
only about 0.25. For k d = 100, however, the range in the 
temperature function increases to about 0.95. The dimen- 
sionless radiative heat fluxes (given in the figure legend) il- 
lustrate how the heat flux decreases as the optical thickness 
is increased. Making n more continuous by increasing the 
number of sublayers increased q r slightly. 


Effect of Total Optical Thickness 

Figure 2 also shows the effect of the total optical thickness 
of the layer. For this figure, the optical thicknesses of each 
of the sublayers are all equal so that for each value of k d , 
k dj = where J is the total number of layers. When k d 
= 0.1, the optical thickness of each sublayer is small enough 



a) Dimensionless coordinate, X=x/D 



Fig. 2 Dimensionless temperature distributions in a composite region 
with the index of refraction increasing from one to three through the 
total thickness. The total optical thicknesses are k d = 0.1, 1, 10, 100, 
and the sublayers have equal thicknesses so that k dj = k d !J : a) number 
of layers, /= 5; b) number of layers, J = 21. 


Effect of Increasing and Decreasing Refractive Indices in a Layer 
For the distributions in Fig. 2, the refractive index increased 
in uniform increments across the entire composite layer. Fig- 
ure 3a shows the effect of a uniform incremental increase 
from n = 1 in the first layer to n = 3 in the central two 



3) Dimensionless coordinate, X=x/D 



b) Dimensionless coordinate, X=x/D 

Fig. 3 Dimensionless temperature distributions in a 10 - layer com- 
posite region with the index of refraction varying in the range from 
one to three. The total optical thicknesses are k d = 0.1, 1, 10, 100, 
and the sublayers have equal thicknesses so that k dj = k d /10: a) 
refractive index increases from one to three and then decreases to one; 
b) refractive index decreases from three to one and then increases to 
three; symbols show solution using exact and <P functions, dashed 
lines are results using diffusion relations for W and <I>. 
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sublayers, and then a uniform incremental decrease to reach 
n = 1 in the final sublayer. The optical thicknesses of the 
sublayers are all equal so that for each sublayer k dj = k d IJ , 
where J = 10. The same total optical thicknesses are used as 
in Fig. 2, k d = 0.1, 1, 10, 100. In Figs. 3a and 3b there is 
the same portion of the total layer thickness at each of the 
refractive indices, as in Fig. 2a. Each of the sublayers in Fig. 
2a has been divided in half and the order of the resulting 10 
layers rearranged. Thus, useful comparisons can be made 
between Figs. 3a and 2a. In Fig. 3a the dimensionless tem- 
perature distributions become flat in the central portion of 
the layer where the refractive index is largest. The profiles 
are asymmetric about the center plane. The effect of increas- 
ing the optical thickness on the temperature distribution through 
the entire laminate is not as pronounced as in Fig. 2. Com- 
paring the heat fluxes in Figs. 3a and 2a shows the rearrange- 
ment of the layers to begin and end with the same index as 
the surroundings has resulted in a small increase in q r . 

Effect of Order of Refractive Indices in Heat Flow Direction 

In Fig. 3b the k d values are the same as in Fig. 3a. The n 
values, however, are interchanged so that the largest refrac- 
tive index, n = 3, is now adjacent to each of the two outer 
boundaries. Then proceeding toward the center, the n de- 
creases in uniform increments to one, such that the n variation 
is symmetric across the entire layer. As compared with Fig. 
3a, the different distribution of n values changes the distri- 
bution of totally reflected energy within the sublayers. There 
is also partial reflection of the incident flux q° rl \ this did not 
occur for the conditions in Fig. 3a where n = 1 in the out- 
ermost sublayers, and therefore, matches the n of the sur- 
roundings. The dimensionless temperature function distri- 
butions have their largest gradients in the center of the layer 
where n has its lowest value. Compared with Fig. 3a, the 
variation of k d has a larger effect. The profiles are again 
asymmetric about the center plane. For small k d , there is a 
significant reduction in the radiative flux being transferred, 
compared to Figs. 2 and 3a. 

Figure 3b also provides a comparison with solutions com- 
puted using the exact functions for <P(X, k d ) and ^(#c D ); the 
exact solutions are plotted for two values of k d using square 
and circular symbols to make them distinct from the dashed 
lines that were obtained using the diffusion functions. The 
comparisons are excellent showing the validity of using the 
very convenient diffusion relations for these basic functions. 

Reciprocal Behavior of Dimensionless Temperatures 

Figure 4 illustrates a general reciprocal behavior for the 
temperature results; an example with five sublayers is used. 
The reciprocal behavior is valid for an arbitrary number of 
layers and for a continuous variation in n. In Fig. 4a, the n 
increases linearly in an incremental fashion from one to three 
progressing from the hot to the cold side. Rather than having 
equal k dj as in Figs. 2 and 3, the k dj have various values to 
illustrate the nature of k dj variations. The physical thicknesses 
of the sublayers is shown as equal, but they could be unequal, 
as the solution depends only on the optical thicknesses k dj 
(this will be discussed in Fig. 4c). In Fig. 4b, both the n values 
and the optical thicknesses of the five layers are reversed in 
their order, relative to the hot side as compared to Fig. 4a. 
The reciprocal result obtained is that rotating Fig. 4a by 180 
deg gives the same dimensionless temperature distributions 
as in Fig. 4b. It was also found that the radiative heat flux 
through the laminate is independent of reversing the order of 
the multiple layers relative to the side with larger incident 
radiation; hence, in Fig. 4b, the heat fluxes remain the same 
as in Fig. 4a. 

Although the radiative behavior depends on the k Dj , rather 
than the physical thicknesses D } of the layers, the temperature 
gradients depend on the Dj. As an illustration, the solution 
in Fig. 4a is given in Fig. 4c, with three of the sublayers having 
small physical thicknesses as might be characteristic of rein- 



a) Dimensionless coordinate, X=x/D 




Fig. 4 Dimensionless temperature distributions in a five - layer com- 
posite region with the index of refraction increasing from one to three 
through the total thickness. The effect is shown of having different 
optical thicknesses for each sublayer. Total optical thickness are k d 
- 1.93 and 19.3: a) refractive index increases from one to three for 
various k dj ; b) variations in both n } and k dj are reversed compared 
with part a; and c) part a with layers having unequal physical thick- 
nesses. 

forcing layers in a ceramic material or a coated part. The 
temperature variation through each sublayer remains the same, 
but in a sublayer with a smaller physical thickness, the local 
gradient of the dimensionless temperature function is in- 
creased. 

Figure 5 shows the effect of having different rates of in- 
crease of the index of refraction within a layer. Consider the 
three solid curves; they are all for a total optical thickness of 
k d = 100. As shown in the legend, the three thicknesses of 
the solid lines correspond to different increases in the refrac- 
tive index within the layer. For the curve drawn with a thin 
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Fig. 5 Effect of the rate of increase of the refractive index within the 
layer for various optical thicknesses. 

line, n — 1.5 throughout the layer. For the curve with a 
medium thickness line, n increases in a stepwise fashion of 
Art, = 0.2 for each sublayer, so the first and final sublayers 
have n = 1.5 and 2.5. For the heaviest line, Arty = 0.5, so 
that n has the range 1. 5-4.0 across the layer. Increasing the 
« variation results in increased curvature of the dimensionless 
temperature function. For k d = 10 there are the same trends, 
this is shown by the three curves drawn with the longer dashes. 
The remaining three curves drawn with the short dashes are 
for a smaller optical thickness, k d = 1. In this instance it was 
found that the curves do not cross and the temperatures in- 
crease as An becomes larger. 

The dimensionless radiative heat flux also has different trends 
as the optical thickness changes from 1 to 100. For k d = 1, 
the radiative flux decreases as the gradient in refractive index 
increases (larger Arty). For an optically thick layer, the be- 
havior is reversed and a greater increase of n with jc increases 
the radiant energy transfer. 

Conclusions 

A method was developed for obtaining temperature dis- 
tributions and radiative heat fluxes in a plane layer with an 
index of refraction and optical density that varies within its 
thickness. The approach that was derived is to analyze a lam- 
inated region with many layers that can each have a different 
refractive index and optical thickness. The interfaces between 
layers are assumed diffuse with the intent of providing infor- 
mation on the limiting case of energy transfer only by radia- 
tion in a ceramic material with layers having various com- 
positions. Radiative energy is incident on each outer boundary. 
The multilayer laminate is in radiative equilibrium (limiting 
case of zero heat conduction) and emits, absorbs, and iso- 
tropically scatters radiation. The method yields results for 
multiple laminated layers with arbitrary indices of refraction 
by using the known dimensionless temperature and heat flux 
functions for a single semitransparent radiating layer with n 
= 1. These functions are combined through the coupling at 
the interfaces to build up solutions for a multilayer laminated 
region. The coupling relations include transmission through 
the interfaces and all internal reflections within the layers. 
The result is a set of simple algebraic expressions for the 
temperature distribution in each of the sublayers and for the 
radiative heat flux through the laminated composite region. 
Sample results were evaluated using the basic functions for a 
single layer with n = 1 obtained by exact numerical integra- 
tion and by diffusion methods. The results using diffusion 
functions were found to be within engineering accuracy. Using 
these functions provides an especially simple way for com- 


puting multilayer results for any distribution of optical thick- 
ness. 

Illustrative dimensionless temperature function distribu- 
tions and radiative fluxes are given for a variety of refractive 
index variations and optical thicknesses of the layers. A rec- 
iprocity relation was found for the temperature function dis- 
tributions. If the order of the indices of refraction of the 
sublayers is reversed as well as their optical thicknesses, the 
temperature distributions become inverted mirror images while 
the radiative heat flux is unchanged. It was found that in- 
creasing the local refractive index in a layer makes the local 
temperature distribution more uniform by means of increased 
internal reflection of energy. The effect that this has on the 
radiative energy transfer depends on the optical thickness of 
the composite. 
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